LEED holography applied to a complex superstructure: 
a direct view of the adatom cluster on SiC(lll)-(3x3) 
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For the example of the SiC(lll)-(3x3) reconstruction we show that a holographic interpretation of discrete Low 
Energy Electron Diffraction (LEED) spot intensities arising from ordered, large unit cell superstructures can give 
direct access to the local geometry of a cluster around an elevated adatom, provided there is only one such 
prominent atom per surface unit cell. By comparing the holographic images obtained from experimental and 
calculated data we illuminate validity, current limits and possible shortcomings of the method. In particular, we 
show that periodic vacancies such as cornerholes may inhibit the correct detection of the atomic positions. By 
contrast, the extra diffraction intensity due to slight substrate reconstructions, as for example buckling, seems 
to have negligible influence on the images. Due to the spatial information depth of the method the stacking of 
the cluster can be imaged down to the fourth layer. Finally, it is demonstrated how this structural knowledge 
of the adcluster geometry can be used to guide the dynamical intensity analysis subsequent to the holographic 
reconstruction and necessary to retrieve the full unit cell structure. 
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I. INTRODUCTION 

The majority of ordered surface structures known 
today has been determined by Low Energy Electron 
Diffraction (LEED)EJ. The strong elastic and inelastic 
interaction of electrons in the energy range 50-500eV in- 
volves a particular sensitivity to the atomic arrangement 
within the outermost layers. In many cases this permits 
a structure determination with a precision of a few hun- 
dredths of an Angstr0m, making LEED one of the pri- 
mary surface crystallography methods. Unfortunately, 
the strong multiple scattering implied by this type of 
interaction also tremendously complicates the theoreti- 
cal analysis of the acquired data. In addition, the real 
space geometry can not be drawn from the intensities 
directly, so that standard quantitative LEED structure 
determinations have to apply a trial-and-error method 
which is frequently supported by structural search pro- 
cedures. Calculated diffraction intensities of a multitude 
of models have to be compared with the experimental 
data until eventually a sufficiently high agreement be- 
tween both is achievedDQ. 

Though more and more advanced experimental and 
theoretical developments have recently given access to 
rather complex surface structureso, it is just this com- 
plexity which at a certain degree inhibits the successful 
application of quantitative LEED. The number of mod- 
els as resulting from the mere combination of all coor- 
dinates of the many atoms in a large unit cell structure 
becomes so huge that it is difficult if not impossible to 
handle. This applies even when for example usingr-au- 
tomated search algorithms in multi-parameter spaceQ as 
the latter has to include the correct model. However, for 
a large unit cell structure our structural imagination is 
frequently unable to even define the type of the correct 
model or the relevant part of the parameter space con- 
taining the real structure, in which a search could then 



be started. Also, methods developed in LEED to deter- 
mine the atomic positions directly still rely on an initial 
good guess of the real structureS. 

The holographic approach represents a revival of the 
hope for the direct disclosure of structural information. 
The idea, ficst developed for the related Photoelectron 
DiffractionEfB, aims at the determination of at least par- 
tial features of the structure when in addition to the 
multiple scattering problem the complexity of the sur- 
face prohibits its full retrieval. In the present case, the 
information provided consists of the local environment 
around an elevated atom in the cell, which might be an 
adsorbate or an intrinsic adatom resulting from surface 
restructuring. Even though only some atomic positions 
at a rather coarse resolution are determined in this man- 
ner, the necessary consistency of the obtained structural 
unit with the complete surface geometry may rule out 
many models directly. Hence, the remaining parameter 
space can be reduced to an extent sufficient to allow the 
application of conventional surface crystallography meth- 
ods. 

The first translation of holographic schemes to-the field 
of LEED as proposed by Saldin and De AndresEI was re- 
stricted to surfaces on which atoms or molecules are ad- 
sorbed in lattice gas disorder. The lacking periodicity 
creates diffraction intensity also outside the sharp sub- 
strate Bragg spots and causes a diffuse intensity distri- 
bution on the screen (for a recent review on Diffuse LEED 
(DLEED) see e.g. ref. [lO| ). This appeared as the natural 
input for the Fourier-like integral transform typical for 
holographic techniques. In the course of subsequent the- 
oretical improvements a proper reconstruction algorithm 
could be established that allowed to circumvent several 
problems complicating the holographic interpretation of 
LEED intensities (see section II). In the present investi- 
gation we use the latest stage of this development which 
allows to construct a reliable image of the complete 3D 
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atomic surrounding of j-the.. elevated atom from data of 
normal incidence aloneETtj. 

However, these theoretical achievements were based on 
the use of diffuse intensity distributions emerging from 
disordered systems while the majority of interesting sur- 
face structures are ordered phases, often with large su- 
perstructure unit cells. Certainly, it would be very ad- 
vantageous to obtain a partial, but direct information by 
holographic means for these onlered phases, too. As we 
briefly demonstrated recentlyEj, the diffraction intensi- 
ties arising from this class of systems may be used as 
input to just the same holographic reconstruction algo- 
rithm as developed for the DLEED case. Two impor- 
tant restrictions for this type of application have to be 
mentioned: there must be only one elevated adatom per 
surface unit cell and the unit cell must have a minimum 
size. The first condition arises from the necessity of a 
unique holographic reference wave as we will outline in 
the third section. The second limit is based on the data 
density available and required. The approximate mini- 
mum size of the-junit cell could recently be estimated as 
a p(2x2) meshli3. An upper size limit is drawn by exper- 
imental factors and the more likely appearance of several 
adatoms per unit cell with increasing unit cell area. A 
(7x7) cell already appears to be too large as discussed in 
section HI. 

Still, a considerable number of ordered reconstructions 
remains open for a holographic analysis. In the present 
paper, our investigations are focused on the first success- 
ful application to an ordered and a priori unknown com- 
plex structure case, the example of the SiC(lll)-(3x3) 
superstructure. This surface phase is of considecable in- 
terest in current crystal growth invesiigationsli3 of the 
pronjising semiconductor material SiCEZI. Previous STM 
workliaEEl had revealed a single large protrusion per sur- 
face unit cell. So, this reconstruction seemed particu- 
larly suited for a first application of holographic LEED 
to ordered surfaces meeting both requirements outlined 
above, i.e. sufRcient unit cell size and the presence of 
a single elevated atom. By comparing the holographic 
images obtained from experimental data and from calcu- 
lated intensities for fictitious models deviating from the 
real surface geometry we illuminate new aspects of the 
validity and possible shortcomings of the new method. 

The paper is organized as follows: in the next section 
we recall the holographic reconstruction algorithm using 
diffuse LEED intensities. Thereafter we describe the re- 
lation between diffuse and discrete intensities and give 
arguments under which circumstances the holographic 
algorithm may readily be applied to conventional spot 
intensities. This is followed by the reconstruction of 
an atomically well resolved image from the experimen- 
tal LEED intensities measured for the SiC(lll)-(3x3) 
phase. Section V shows that the spatial depth accessible 
by the method is rather large, which allows to deter- 
mine such important features as the stacking sequence of 
deeper layers. In section VI, we illuminate the role of pe- 
riodic vacancies within the unit cell acting as additional 




FIG. 1. Schematic display of the holographic interpretation 
of the adatom scattering: a) electrons finally scattering at the 
adatom form the reference wave R, those subsequently hit- 
ting one substrate atom represent the kinematic object wave 
O (solid lines). The dashed lines display possible multiple 
scattering events providing dynamic contributions to the ref- 
erence and object wave (see text for details), b) Pronounced 
forward scattering at the beam-splitter indicated by difi'erent 
length of the arrows in different directions. 

holographic reference waves. Then we address the issue 
of intensities arising from substrate relaxations such as 
buckling and finally discuss the use of the holographic in- 
formation for a complete surface structure analysis in the 
case of SiC(lll)-(3x3), whose precise realjSpace struc- 
ture is described in more detail elsewhereEil'cSl. 



II. THE HOLOGRAPHIC RECONSTRUCTION 
ALGORITHM 

The holographic approach in DLEED makes use of 
the fact, that all measurable diffuse intensity outside the 
sharp substrate Bragg spots necessarily has been caused 
by at least one scattering event at one of the disordered 
adsorbates on top of the (unreconstructed) crystal: scat- 
tering exclusively within an ideal bulk-terminated sub- 
strate can only lead to diffraction intensities at Bragg 
spot positions. In that sense the adsorbate atom can be 
viewed as a prominent scatterer which, acting as a beam- 
splitter, provides a natural separation of all scattering 
paths as depicted by the solid lines in Fig. 1(a): electrons 
whose final scattering is by an adsorbate form the refer- 
ence wave i?(k), while those scattered subsequently by 
substrate atoms before reaching the detector provide the 
object wave 0(k)l£l (where k = (ky, k±) is the wavevector 
of the detected electron, with the components k|| parallel 
to, and fcj_ perpendicular to the surface). This allows to 
interprete the diffuse intensity as the interference pattern 
of these two contributions. Hence, the local surround- 
ing of the beam-splitter atom should bcj^tractable by a 
phased 2D Fourier transform of the datafl'Ll. However, the 
above interpretation does not include that multiple scat- 
tering adds unwanted contributions to both reference and 
object wave as indicated in Fig. 1(a) by the dashed lines. 
A considerable improvement in the image quality could 
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be achieved by combining several DlLEED patterns mea- 
sured at different electron energiesES. The corresponding 
multi-energy reconstruction algorithms include a 3D in- 
tegral transform and try to single out the contributions 
due to the kinematic object wave, suppressing the un- 
wanted effects caused by the multiple scattering of the 
low energy electrons. 

Yet, the pronounced forward scattering of the beam- 
splitter led to only a selective appearance of the atoms in 
the reconstructed local adsorption geometries, depending 
on whcthcr|-they were located within the forward scat- 
tering conecj, cf. Fig. 1(b). The implied necessity of 
combining several (at least two) data sets taken at dif- 
ferent angles of incidence to deduce the complete 3D sur- 
rounding of the beam-splitteiEj could be overcome with 
the introduction of an improved. reconstruction algorithm 
proposed by Saldin and ChenHl. 

This Compensated Object and Reference wave Recon- 
struction by pj Energy- dependent Cartesian Transform 
(CORRECT)EI allows the calculation of the real space 
distribution around the adsorbate \B{r)\'^ (where r = 
(r||,z) is a position vector relative to the origin at the 
adsorbate with components ry parallel to, and z perpen- 
dicular to the surface) via the following expression: 

B{r)= ff \f if(k||,A;^;r)x(k||,A:^)e-(^'='-'=-^)&x 

(2.1) 



reference wave. In a zeroth order approximation it can 
be written 



Note, that in contrast to previous reconstruction al- 
gorithms that performed the involved 3D integral in a 
polar coordinate system (angle and energy), the data in- 
put is provided on a cartesian grid (k|| ,k±), which will be 
of importance when discussing the step towards ordered 
superstructures in the next section. 

The transform does not operate directly on the mea- 
sured intensities H, but rather on a contrast-enhancing 
and normalizing function 



with 



g(k||,fcx)-ga.(k||) 
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(2.2) 



(2.3) 



It has been shown theoretically that the use of such a 
X-function helps to partially remove the self-interference 
terms |i?(k)|^ and |0(k)|^ in the DLEED intensity, which 
give rise to spurious high values of the reaLspace distribu- 
tion |i?(r)|^ in the vicinity of the origintJ. Additionally, 
X has been designed in such a way as to suppress mod- 
ulations in the DLEED patterns that—arise from some 
partial ordering among the adsorbatesEj. 

The last part in the expression to be described is the 
integral kernel which corrects for the anisotropy of the 



Km,kr,r) = 



/a(k. ■r) + C 



(2.4) 



Here /a(ki • f) is the atomic scattering factor of the ad- 
sorbate, ki the direction of electron incidence, and C 
the so called kernel constant (which we take to be real), 
and which represents an isotropic approximation to the 
backscattering by the substrate prior to scattering by 
the adsorbate. Optimizing the value of C provides ac- 
cess to those atoms of the local adsorption geometry that 
lie outside the forward scattering direction of the beam- 
splitteic3. This allows the retrieval of the complete 3D 
surrounding of the latter from data of normal incidence 
alone. The algorithm in this present forn\iias.been shown 
to give reliable images usipg[-tiieoreticaLJ'E3, as well as 
experimental DLEED datalll'tj. 



III. SPOT INTENSITIES VERSUS DIFFUSE 
DISTRIBUTIONS 



The original holographic reasoninga was based on the 
assumption that only one beam-splitting adsorbate atom 
is present on the substrate surface. With several such 
adsorbates, each time in the same local structure but 
without long range order among them (lattice gas dis- 
order), intensities simply add up in the low coverage 
limit leaving the resulting diffuse distribution practically 
unchangeal3'E3. A different situation emerges with the 
onset of order at higher coverages and/or upon thermal 
annealing: additional modulations in the DLEED pat- 
tern are created that eveiittjially cause the breakdown of 
the holographic algorithmll^. 

For the case of a completely ordered superstructure 
of such adsorbates, the modulations caused by the lat- 
tice factor concentrate the diffuse intensities to a series 
of discrete superstructure - or fractional order - spots 
when the unit mesh of the adsorbate layer is larger than 
that of the crystalline substrate. However, the simulta- 
neous extinction of diffuse intensity between the spots 
as caused by destructive interference between the waves 
originating from different adsorbate-substrate clusters, 
does not remove the crystallographic information wanted: 
the energy dependence of the superstructure spot inten- 
sities is the same as the one displayed by the correspond- 
ing k|| -positions in a diffuse distribution resulting from 
a disordeced adlayer in the equivalent local adsorption 
geometrycS. The only restriction is that scattering be- 
tween such clusters has to be negligible, a condition sat- 
isfied even for relativeliz-small superstructures when using 
normal incidence dataou. 

So, even though the perfect order among the adsor- 
bates significantly reduces the amount of available data, 
the few remaining intensities are not masked by disturb- 
ing modulations as in the case of partial disorder inhibit- 
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ing the final ky -integration in equation (2.1). The su- 
perstructure spots can be thought of as sampling the 
DLEED intensity distribution of the corresponding lat- 
tice gas on a finite grid. Therefore, as suggested earlierE3, 
a DLEED holographic algorithm may in principle be ap- 
plied to such ordered superstructure systems with the 
only difference of a reduced density of input data in k|| . 
This makes more apparent, why the CORRECT algo- 
rithm is so particularly suited for the extension to ordered 
phases: the data is provided on the appropriate cartesian 
grid and only normal incidence is required. Interestingly, 
earlier investigations on the information content of dif- 
fuse intensitiesEa, as well as on the minimum data base of 
the algorithmEa, showed, that the continuous diffraction 
distribution resulting from disordered atomic adsorbates 
is already sufficiently described when using a (3x3) sam- 
pling grid. Information on a denser grid is largely redun- 
dant. Hence, there are no drastic changes to be expected 
when making the transition from disordered systems to 
phases with large superstructure cells like the (3x3) re- 
construction of SiC(lll) of the present paper. This al- 
lows us to apply the algorithm developed for DLEED 
without any modifications. However, it should be em- 
phasized, that a further reduced data base in connection 
with superstructures smaller than a (2x2) can lead to 
aliasing effects iH-|the Fourier-like transform due to insuf- 
ficient samplingtj. 

The application of LEED holography to ordered sur- 
faces involves several practical advantages. For the 
diffraction process it is irrelevant, whether the beam- 
splitter is an externally adsorbed atom or intrinsically 
belongs to the surface. Thus, besides ordered adsorp- 
tion systems now also ordered substrate reconstructions 
can be investigated. Additionally, the measurement of 
discrete spot intensities is much less delicate than that 
of the diffuse intensities which are comparatively weak. 
The high signal-to-noise ratio of the bright spots allows 
easy subtraction of contributions due to thermal dif- 
fuse scattering. Also, at higher energies fractional spot 
intensities are not that much influenced by cross-talk 
from the Iwgift substrate spots as is the case for diffuse 
intensitiesEj'ElJ. Furthermore, holographic LEED seems 
also suitable to tackle larger unit cell reconstructions: 
the high number of fractional order spots generated in 
these cases provides a fine sampling grid and ensures the 
proper working of the integral transform. However, prac- 
tical reasons also commend an upper limit for the unit 
cell size as the increasing number of closely spaced spots 
impairs a proper data acquisition, especially at higher 
energies where more and more spots appear and weak 
spots are disturbed by their brioht .neighbours. A unit 
cell such as the (7x7) on Si(lll)E]L3 is probably already 
too large from an experimental point of view as the ac- 
cessible energy ranges become too small. 

In addition, it becomes more and more unlikely that 
such a large unit cell contains only one elevated adatom 
(the Si(lll)-(7x7) actually contains 12 adatoms). This 
would violate the strongest restriction of the technique 



at its current stage, i.e. the condition that only a single 
beam-splitter is allowed within each unit cell. Several 
such prominent atoms per unit cell would lead to inter- 
mixing of their respective contributions as will be demon- 
strated further below. This is all the more problematic, 
since the actual number of elevated adatoms is just one 
of the quantities sought in the structure analysis of an 
a priori unknown surface (even though STM might help 
as in the present case). Future efforts in methodologic 
improvements should hence be directed to overcome the 
multiple beam-splitter problem, which did not occur in 
the previous applications to simpler diffuse or ordered 
systems. As a consequence, until there is a proper theo- 
retical description of the detailed influences on the recon- 
structed images, the systems to which holographic LEED 
is to be applied have to be chosen with considerable care. 



IV. RECONSTRUCTION USING 
EXPERIMENTAL DATA 

SiC is a material that displays most suitable electronic 
properties which have made it a promising candidate for 
high power and high frequency devices. Particularly, 
the (3x3) phase of the SiC(lll) surface has drawn con- 
siderable interest in the last years ,|-caused by the ob- 
served crystal growth improvementE3, that is achieved 
when this paconstruction is stabilized under highly Si-rich 
conditionals. Its complexity, which can ahead v_be_jia- 
duced from an extensive debate in the literatureErEElEj, 
had hitherto prevented a detailed structure analysis us- 
ing trial-and-error methods. However, the high number 
of fractional order spots caused by such a comparatively 
large surface unit cell makes this phase an ideal candidate 
for a holographic investigation in view of the reasoning 
outlined above. 

The cubic 3C-SiC polytype was chosen, since its 
(111) oriented surface exposes only one definite stack- 
ing sequencetj, i.e. there is no coexistence of domains of 
different orientation, which would have to be expected in 
the case of hexagonal polytyaes with different layer stack- 
ings possible at the surfaceEj, and which would certainly 
complicate if not inhibit the interpretation of the recon- 
structed images. Additionally, there is strong evidence 
from comparison of experimental LEED intensitiesEj, as 
well as from DFT test calculationsEJ, that the atomic 
structure of the (3x3) surface phase itself is rather in- 
dependent of the sample polytype. So, results obtained 
for 3C-SiC(lll) can be expected to hold also for other 
polytypes. 

LEED I(V)-curves of the sharp diffraction pattern were 
measured in the energy range 50-300 eV using normal 
electron incidence. Details on the data acquiaitioji and 
sample preparation will be published elsewherecj'EJ. The 
low diffuse background and noise level allowed the record- 
ing of 14 fractional order beams closest to specular re- 
flection, which are symmetry-inequivalent at normal in- 
cidence. Providing the measured intensities as input to 
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FIG. 2. Recovered local geometry of the SiC(lll)-(3x3) 
structure using experimental data in the energy range 
50-300 eV and kernel constant C = 2.7 A. The maximum 
noise level in the image is 48 % of the maxima denoting the 
atom positions (noise cut-off: 25 %). For details on the dis- 
play procedure, see Section IV. The inset displays a schematic 
of the retrieved adcluster geometry including chemical bonds 
and the approximate layer distances as determined by holo- 
graphic LEED. 



expression (2.1) resulted in the 3D image displayed in 
Fig. 2: the real-space distribution |i3(r)|^ is calculated 
on a grid of 0.2 A resolution inside a cylinder of depth 
6.0 A and a lateral radius 3.0 A, which is consistent with 
estimates on the lateral validity of the algorithmtj. Small 
spheres are drawn at the grid points, indicating the re- 
constructed real-space intensity by their diameter which 
scales linearly with the intensitv—As. ipointed out in pre- 
vious holographic investigationslliriljO, this type of dis- 
play permits a quick understanding of the essential fea- 
tures of the structural unit determined holographically 
and will therefore be used in all figures included in the 
present paper. 

The origin of the coordinate system is defined by the 
beam-splitter, which is artificially added in the image as 
a black sphere to facilitate understanding. The highly Si- 
rich conditions under which the (3x3) phase is observed 
suggest this beam-splitter to be Si, the scattering factor 
of which is cons equ ently used for the computation of the 
integral kernel (2.4). However, the zeroth order approxi- 
mation of the latter is most sensitive only to the essential 
form of the atomic scattering factor, which is very sim- 
ilar for most elements. Using carbon as a beam-splitter 
in the computation consequently did not change the re- 
sulting images considerably. The kernel constant C in 
this expression is optimized such that all atoms in the 
geometry appear with approximately equal brightnesses. 



The highest disturbing intensity at non-atomic positions 
(henceforth referred to as noise level) is with 48 % of the 
overall maximum value at an unprecedented low level. 
This has to be attributed to the - in comparison to the 
BLEED case - much better quality of conventional LEED 
I(V) data and the increased energy range available. 

The image allows the unambiguous identification of 
the local adcluster geometry formed by an adatom sup- 
porting trimer and two further atoms vertically below 
the beam-splitter (see inset in Fig. 2). The rough layer 
distances of 1.3 A (adatom-trimer), 1.3 A (trimer and 
first lower atom) and 2.0 A (between lower atoms) corre- 
spond surprisingly well with the (7x7) DAS model of the 
Si(lll) surfaca23. This already indicates that probably 
the complete retrieved geometry corresponds to Si atoms 
on top of the SiC substrate. Note, that the distorted form 
of the trimer atoms is an effect of the scattering factor in 
connection with the zeroth order approximation of this 



property in the intcgr al kernel (|2j)El. However, neither 
the obtained spatial resolution, nor the exact position of 
the atoms inside the geometry are the primary object or 
strength of the holographic analysis: it is rather the di- 
rect and quick idea of a structural unit belonging to the 
investigated surface. 

The obtained tetramer formed by the adatom and the 
supporting trimer is typical for hexagonal semiconduc- 
tor surfaces. Its unambiguous determination in the holo- 
graphic image proves that only one of the two possible 
orientations rotated by 60° with respect to each other 
is present on the surface. This already excludes do- 
mains of differently rotated, i.e. coordinated clusters, 
as can also directly be deduced from the pronounced 
threefold pmrnetry of the measured fractional order spot 



intcnsiticstil[ill.|-4t further rules out the model proposed 
first by KaplanEJ of a (3><r3j) mesh which in close analogy 
to the (7x7) DAS modelEa contains two such tetramers 
per surface unit cell, which in turn would necessarily be 
differently oriented. It should be noted, that this model 
was also inconsistent with STM investigations, which 
clearly revealed only one elevated protrusion per unit 
cell, thus-sttpngiy favouring models including a single 
tetramerllac2l. Since the atomic beam-splitter has to be 
identified with the top adatom of this tetramer, exactly 
these results ultimately enabled the application of LEED 
holography to this structure: the obligatory uniqueness 
of the beam-splitter excludes DAS-like models with two 
tetramers per surface unit cell from the class of systems 
accessible under the current state of theory. 

We should recall now that we are dealing with an a 'pri- 
ori unknown structure. Although the low noise level in 
the image may appear very convincing, it has to be recog- 
nized that the strong multiple scattering combined with 
the anisotropic scattering factors for low energy electrons 
may lead to serious artefacts in the images that would 
not easily be distinguishable from real atoms. Since the 
multi-energy algorithms developed for holographic LEED 
can only suppress, but not completely eliminate these ef- 
fects, it is often advisable to vary the used input energy 
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range. In view of the fact that scattering factor anomahes 
may sometimes even lead to an increased image quahty 
when reducing the number of included energies, the sta- 
bility of the obtained result under such variations can 
significantly increase the confidence in the deduced lo- 
cal geometry. Dividing the experimental data into vari- 
ous subsets always resulted in equivalent images, which 
strongly confirms the structural unit determined. In gen- 
eral, one has to admit that the weak scattering power of 
light elements like Si and C might provide a favourable 
case for holography as multiple scattering contributions 
are expected to be smaller than for example for transition 
metal crystals. However, recent results for the system 
O/Ni(001)-p(2x2)y make us believe, that the developed 
algorithm does also work for stronger scattering materi- 
als. Now, although the experimental result appears con- 
vincing a test of the validity and sensitivity of the method 
seems adequate, which is presented in the next sections 
using simulated intensities from various fictitious models. 

V. VERTICAL SENSITIVITY AND STACKING 
OF DEEPER LAYERS 

The simplest model consistent with the atomic posi- 
tions obtained from the holographic reconstruction would 
be just an adatom on top of a SiC bilayer. Yet, such a 
model appears improbable since it would not account for 
the strong silicon enrichment at the surface as detected 
by earlier Electron Energy Loss (EELS) and Auger Elec- 
tron Spectroscopy (AES) results^, which we discuss in 
detail elsewhereEHl. Assuming therefore the tetramer as 
the essential structural unit - presumably formed by Si 
atoms in view of the EELS and AES results - the first 
question to be verified is its position on the underlying 
substrate. The simplest possible solution would be to 
directly place it somehow on top of the SiC sample, to 
which the two further atoms showing up in the holo- 
graphic reconstruction would then belong. Given the 
bilayer stacking sequence of this material in the [111] 
direction it is, however, impossible to find any location 
consistent with two atoms directly on top of each other as 
predicted by Fig. 2. Even though the holographic method 
is most reliable in just this direction vertically below the 
beam-splitter, in which atoms show up already without 
the scat tering factor compensation by the integral ker- 
nel there is yet no definite certainty on the limit 
of the algorithm's validity for deeper layers: all previous 
investigations with BLEED data had dealt with rather 
simple test structures, which were already completely de- 
termined by the atomic positions in the first two layers. 
Hence, the now performed calculation of |i3(r)|^ up to a 
depth of 6.0 A raises concerns whether the lowest lying 
atom identified at 4.6 A might already be outside such a 
limit. 

Consequently, for the moment we focus on only the 
tetramer and the third layer atom of the holographic 
image. This arrangement suggests a cluster position in 




FIG. 3. Same as Fig. 2, but using simulated I(V) curves of 
a simplified Li/Tsong model as described in section V. The 
electron energy range was 146-300 eV, the kernel constant 
C = 5.0 A and the maximum noise level is at 46 % of the 
maxima at the atom positions (noise cut-ofF: 25 %). The in- 
set displays the atomic arrangement in the model assumed to 
calculate the intensities used for the holographic reconstruc- 
tion. 

which each of the three trimer atoms is located in a hol- 
low site of the topmost hexagonal bilayer of the SiC sub- 
strate. Depending on whether this position is fourfold 
coordinated, i.e. on top of a carbon atom in the sub- 
strate bilayer (T4 site) or truely threefold coordinated 
(H3 site), the tetramer would have to be oriented as in 
the SiC substrate bilayers or rotated by 60° , respectively. 
The adatom would then reside on top of a Si atom in the 
bilayer, which corresponds to the third layer atom of the 
holographic image. No further atom would be present 
2 A below the latter in either site geometry. Note, that 
this geometry corresponds to the model proposed by Li 
and Tsong on the basis of their STM workEfl, which as- 
sumed such coordinated tetramers in a (3x3) periodicity 
directly on top of the SiC. In order to verify whether 
the atom additionally appearing in Fig. 2 is an algorith- 
mic artefact or not, we thus simulated theoretical LEED 
I(V) curves of an adcluster model for the identical num- 
ber of fractional order beams and the same energy range 
as in the experiment-Xdetails of these calculations will be 
published elsewhercEHl) . In order to focus on the depth 
information available from LEED holography we chose 
the fourfold coordinated trimer atom positions, yet arti- 
ficially expanded the distance between adcluster and sub- 
strate to push the third layer atom to a position 3.2 A 
below the beam-splitter. A schematic view of this ge- 
ometry can be seen in the inset of Fig. 3. The resulting 
holographic image is displayed in Fig. 3, showing the ex- 
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pected tetramer unit plus the third layer atom, but also 
consistently not indicating any sign of possible artefacts 
vertically below the adatom. Only the carbon atoms of 
the substrate SiC bilayer are still within the reconstruc- 
tion volume, but do not appear in the reconstructed im- 
age (Fig. 3). However, one has to consider that in the 
geometry chosen they are 4.3 A apart from the beam- 
splitter and 1.9 A off the vertical axis, and in addition 
represent comparatively weak scatterers. This distance 
- not in forward scattering direction - probably repre- 
sents the detection limit at least for a weak scatterer. 
However, in view of the absence of artefacts in the holo- 
graphic image obtained for our test model, the presence 
of two atoms vertically below the adatom as indicated 
in the real space reconstruction from the experimental 
data has to be assumed correct. Furthermore, the pro- 
nounced appearance of the lowest atom might indicate 
that it is silicon since a weakly scattering C atom should 
not be detectable at that depth. Our test case resulting 
in the image shown in Fig. 3 thus clearly rules out the 
possibility of the Li/Tsong model, whose bilayer stacking 
sequence results in fourth layer atoms off the vertical axis 
and which is thus incompatible with the lowest atom in 
the image obtained from the experimental data. 

This further emphasizes the importance of the retrieval 
of the two deepest atoms in the local geometry. Since 
only small deviations from the bulk positions are usually 
to be expected in such deep layers, the location of each of 
these atoms uniquely determines the complete stacking 
sequence of the corresponding entire layer. Hence, even 
though the obtained structural unit itself may contain 
only a small number of atoms, its consistent embedding 
into the surface unit cell subsequently can reveal a quite 
important further fraction of the investigated surface. 

VI. CAN A VACANCY ACT AS A 
BEAM-SPLITTER? 

The confirmation of the lowest atom inside the revealed 
structural unit, whose on top stacking is inconsistent with 
a SiC bilayer at the very surface, necessitates to include 
an additional Si adlayer in the crystallographic model. 
Such an adlayer between tetramer and substrate had al- 
ready been included in the original DAS-like model, since 
the EELS and AES results indicated the strcpg presence 
of Si-Si bonding in the highly Si-rich surfaceEj. In order 
to bring this otherwise very reasonable model in accor- 
dance with the STM data described above, Kulakov et al. 
proposed the_a.bsence of one of the two tetramers per sur- 
face unit celllla. In the language of the DAS-type models 
the top atom of the trimer represents an adatom on top of 
a Si bilayer. In the DAS-unit cell one of the two adatoms 
is located on a piece of bilayer in faulted orientation as 
indicated in Fig. 4(a). One would expect an energetical 
difference between the adcluster which follows the sub- 
strate stacking direction and the one, which introduces 




< (3x3)-unit cell > 



FIG. 4. Schematic side view of different models for the 
SiC(lll)-(3x3) reconstruction in a projection parallel to the 
[110]-plane (Si atoms are depicted by large spheres, C atoms 
by small darker spheres. Bonds within the projection plane 
are drawn as single lines, double lines represent two bonds 
pointing out of the projection plane by +60° and -60°, re- 
spectively.) a) DAS-model containing two adatoms and one 
cornerhole per unit cell as proposed by KaplarO. b) Single 
adatom model containing two cornerholes with a local stack- 
ing fault in the Si-bilayer fragment underneath the adatom; 
deri\jed from the model by Kulakov et al. without the stacking 
faultu, c) Single adatom-trimer-cluster residing on a complete 
monolayer in (1x1) periodicity with all cornerholes filled. The 
bilayer fragment underneath the adatom again represents a 
local stacking fault. 

a local stacking fault in the adlayer. Thus it is plausi- 
ble that in the end exclusively the more favourable type 
would be present on the surface. Such a model, hence 
including only a single tetramer with definite orientation 
in each (3x3) unit cell, could explain not only the sin- 
gle protrusion in the STM images, but also the threefold 
symmetry of the LEED pattern. What remains, is the 
question which of both orientations is actually realized. 
Note, that the orientation of the adcluster can be deduced 
from a comparison with a previousJjiEEIj) analysis of the 
(1x1) phase on the same sampleEilSH However, we 
demonstrate a verification of this assignment using test 
calculations, a method that could generally be applied in 
cases where no independent analysis of the substrate is 
available. 

Both Kulakov-type models, with the adcluster either 
introducing or not introducing a local stacking fault, 
comply with the holographic image from the experimen- 
tal data, when identifying the upper of the bottom two 
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FIG. 5. Same as Fig. 2, but using simulated I(V) curves of 
a simplified Kulakov model as described in section VI. The 
electron energy range was 146-300 eV and the kernel con- 
stant C = 5.0 A. a) geometry including cornerholes in the Si 
adlayer, cf. Fig. 4(b): maximum noise level at 76 %, b) ge- 
ometry with filled cornerholes in the Si adlayer, cf. Fig. 4(c): 
maximum noise level is at 29 % (noise cut-off: 25 %). 

atoms as belonging to the Si adlayer and the other one 
as a Si atom of the substrate's topmost bilayer. Yet, 
the subsequent interplay between LEED holography and 
conventional LEED can do better than that: when recon- 
structing images using theoretically simulated data, the 
orientation of the substrate inside the given coordinate 
system is known. Depending on the resulting orienta- 
tion of the adcluster in the image - or to be exact, the 
tetramer of four atoms representing the Si bilayer un- 
derneath the adatom - when choosing one of the two 
equivalent beam assignments in the LEED pattern, its 
orientation with respect to the bulk can be deduced even 
without seeing the latter in the reconstructed image it- 
self. From the result shown in Fig. 3, we therefore know 
the orientation of an unrotated tetramer. Comparing this 
with the holographic image obtained from the experimen- 
tal data (cf. Fig. 2) we find that the adcluster geometry 
involves a local stacking fault of the Si bilayer as shown 
in Fig. 4(b). Hence, of all the previously existing models 
of the SiC(lll)-(3x3) phase, LEED holography would 
only be fully consistent with the Kulakov model in the 
local stacking fault version. 

In order to further verify this conclusion, we simulated 
LEED I(V) curves for this model, whereby, however, 
small deviations from bulk-like positions as for example 
induced by dimerization were not considered. Surpris- 
ingly, the corresponding holographic image displayed in 
Fig. 5(a) is of considerably worse quality than the previ- 
ous results. Although all five atoms of the expected struc- 
tural unit show up, their overall configuration is badly 
distorted and high noise in form of three concentrated 
artefacts prevents the unambiguous distinction between 
real atoms and false contributions. Since the image from 
ideal theoretical data should only be better than the one 
from the experiment, the situation in the latter has some- 



how to be more favourable for LEED holography than so 
far assumed, which, of course, demands clarification. 

As mentioned above, a severe source of an algorithmic 
breakdown at its current level of development is given 
by the multiple beam-splitter problem. Therefore, we re- 
considered the atomic arrangement of the Kulakov model 
under this point of view. Since it was derived in close 
analogy to the DAS model, its Si adlayer is not com- 
pletely closed but contains vacap;ies to relax the stress 
induced by the lattice mismatchE£l. These so called cofc. 
nerholes appear just like in the Si(lll)-(7x7) structured 
and break the (1x1) periodicity as much as the identi- 
fied beam-splitting atom on top of the tetramer. Con- 
sequently one might also concede them the same holo- 
graphic interpretation. In the sketch of the Kulakov 
model in Fig. 4(b) it can be seen that the vacancy is 
even enlarged by the removal of one adcluster from the 
DAS model, cf. Fig. 4(a). It may appear difficult at first 
glance to imagine a vacancy as a possible beam-splitter, 
but speaking in terms of missing wave contributions to 
achieve the destructive interference corresponding to a 
perfect (1x1) adlayer helps to understand its influence 
on the fractional order beams. This would be equivalent 
to replacing the vacancy by a pseudo-adatom with the 
same dynamic scattering behaviour. 

From this point of view, we would have to conclude 
that the Kulakov model contains various distinct beam- 
splitters per surface unit cell, whose respective wave con- 
tributions could consequently interfere and completely 
prohibit a holographical interpretation. However, the 
strong damping of the low energy electrons makes us 
hope, that the dominant contribution in the recon- 
structed image is due to the most elevated, periodicity 
breaking atom in the surface unit cell and that the ex- 
istence of further extra atoms or vacancies in the super- 
structure unit cell leads "only" to image disturbances al- 
though they might be considerable. This interpretation 
would help to understand why Fig. 5(a) basically shows 
the local environment of the top tetramer atom plus arte- 
facts and some distortions that may then be due to the 
cornerhole vacancies. To test this line of thought, we 
simulated LEED I(V) data of exactly the same model as 
before, but filling the vacancies with Si atoms at bulk- 
like positions in the Si adlayer as shown in Fig. 4(c). The 
resulting image in Fig. 5(b) is of impressive clarity and 
contains all essential features of the result obtained with 
the experimental data. We take this as a strong indi- 
cation of the correctness of our reasoning, although we 
want to stress again that there is yet no proper theoret- 
ical treatment of the multiple beam-splitter difficulty in 
LEED holography. It should further be noted that in 
our test model, cf. Fig. 4(c) we only filled the Si mono- 
layer. The adatom supporting trimer atoms are still not 
repeated with the bulk periodicity and thus break the 
(1x1) periodicity, too. 
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VII. INFLUENCE OF SUBSTRATE 
RECONSTRUCTIONS 



What had started as a pure necessity to ensure the cor- 
rect working of the holographic algorithm, subsequently 
turned out to be the last required piece for the solution of 
the (3x3) puzzle. Since LEED holography seemed only 
fully consistent with a Kulakov-derived model contain- 
ing one tetramcr in local stacking fault orientation on 
a closed Si adlayer without cornerholes, a careful recon- 
sideration then showed that indeed there had been no 
other reason for including the cornerholes at first hand 
but the sole analogy to the DAS model. The situation for 
SiC(lll)-(3x3), where the Si adlayer shows an intrinsic 
lattice mismatch of 20 % with respect to the underlying 
substrate, might however require a different form of relax- 
ing the lattice strain under simultaneous dangling bond 
saturation than the corncrhole and dimerization prin- 
ciple underlying the honioepitactic Si(lll)-(7x7). As 
a furtiLML.hint, the obtained STM images of the (3x3) 
phaseEja never showed comparably stroBg cornerhole 
depressions as visible on the silicon surfaceEj. 

Consequently, the thus most favoured model with filled 
cornerholes was input to a refining LEED and Density 
Functional Theory (DFT) analysis. Even though the 
holographic results had considerably reduced the multi- 
parameter space for the trial-and-error search of both 
methods, it should however be emphasized that the re- 
maining structure determination was still a tremendous 
task: there is a qualitative difference between a coarse 
local beam-splitter surrounding that depicts a small frac- 
tion of the huge surface unit cell and a detailed variation 
of all involved atomic positions on a dense grid in steps 
of a few hundredths of Angstr0ms. It was in this respect 
most gratifying that both analyses independently yielded 
the same full (3x3) structure by input of the holograph, 
ically recovered cluster: the resulting final twist modelEll 
can indeed essentially be described as a SiC substrate 
with a strongly buckled Si adlayer without any corner- 
holes plus one tetramer per surface unit cell consisting of 
a trinier and one adatom. These trimer atoms and the 
Si adlayer below locally resemble a Si bilayer in stack- 
ing fault orientation. A more detailed description of the 
exact made.l .and both analyses involved will be given 
elsewhereCJO. 

Yet, we have to realize that the bond optimizing re- 
laxations inside the adlayer also contribute to the super- 
structure spot intensities and might consequently affect 
the working of the holographic reconstruction algorithm: 
each atom, that has left its (Ixl)-like position, has in 
principle to be regarded as a possible additional beam- 
splitter in view of the discussion in section VI. Therefore, 
as a final test, we used the I(V) curves calculated for the 
exact geometry of the optimized model as input to the 
holographic algorithm. The resulting image is displayed 
in Fig. 6 and shows exactly the same local adatom en- 
vironment as the experimental data, which we - even in 




FIG. 6. Same as Fig. 2, but using simulated I(V) curves of 
the final twist model as described in section VII. The electron 
energy range was 65-300 eV, the kernel constant C = 0.75 A 
and the maximum noise level is at 51 % (noise cut-off: 25 %). 

the presence of buckling - take as a final proof of the 
validity of our new method. In accordance with the ex- 
perience made recently with other, simpler structureslij, 
the effect of the slight deviations from bulk-like positions 
on the reconstructed image is apparently negligible and 
leads only to an increased overall noise level, which can 
be seen comparing Fig. 5(b) and Fig. 6, whose under- 
lying structure differs exclusively by just these substrate 
relaxations. All this becomes more understandable, when 
recalling, that it is again only the difference in the outgo- 
ing wavelets arising from buckled atoms that can act as 
a conduit for diffraction intensity in the fractional order 
beams. The (additionally damped) contributions due to 
these shifts can hence be regarded small with respect to 
the major rupture of periodicity caused by the introduc- 
tion of a completely new and elevated atom such as the 
adatom in the present structure. In this context, it is also 
important to notice, that the majority of such shifts is 
far below the resolution capabilities of LEED holography 
at its present stage, which can typically be stated to be 
« 0.6 A. 



VIII. CONCLUSION 

In the present paper we described in detail the contri- 
bution that holographic LEED can provide, when apply- 
ing it to a complex superstructure. Using a holographic 
interpretation of fractional order spot intensities, a 3D 
image of the local geometry around an elevated, periodic- 
ity breaking adatom can be retrieved. The structural unit 
thus obtained has to be consistent with the real space 
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atomic structure and can be used to considerably reduce 
the multi-parameter space and possibly enable the trial- 
and-error search of geometry optimizing methods like 
quantitative LEED and DFT energy minimization. We 
exemplified this for the case of the SiC(lll)-(3x3) phase, 
where the holographically derived adcluster directly rules 
out the majority of all previously existing models of this 
surface and whose geometry has now been fully confirmed 
by the final twist model obtained independently by con- 
ventional LEED and DFT. This application additionally 
marks the first example that a holographic inversion of 
LEED data actually played a crucial part in the deter- 
mination of a complex and a priori unknown structure. 

We have illuminated the power and the limits of 
this new holographic LEED method for ordered sur- 
faces. Even in the best of all cases the obtained im- 
age is restricted to the local geometry around the promi- 
nent beam-splitter. Only for very simple surfaces, this 
uniquely determines all atomic positions inside the unit 
cell. Furthermore, exclusively for such simple test sur- 
faces has this still developing method been thoroughly 
tested so far. There, some severe problems like the mul- 
tiple beam-splitter problem encountered in the present 
investigation do usually not arise and have therefore not 
yet been theoretically treated. Consequently, the sys- 
tems, to which holographic LEED is to be applied, have 
for the moment to be chosen with great care. 

Nevertheless, regarding the immense problems that 
quantitative LEED and DFT face with complex, large 
superstructure systems, every directly obtainable pre- 
information is highly welcome. In this view, the possibil- 
ity of a holographic approach raises hopes that a new ally 
in the structure analysis of these ordered phases has been 
found. It is particularly its direct simplicity that already 
now renders holographic LEED such an ideal supplement 
to its established trial-and-error brethren, even though 
the young method has still a long way to go. 
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